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A general consideration is presented for the magnetic interaction at an interface between a per- 
ovskite manganite and other transition metal oxides. The latter is specified by the electron number 
n in the d3,^2_^2 level as {d-j,^2_,.2)". Based on the molecular orbitals formed at the interface and 
the generalized Hund's rule, the sign of the magnetic interaction is rather uniquely determined. The 
exception is when the d^^2_^2 orbital is stabilized in the interfacial manganite layer neighboring to a 
{dsz^-r^y or (d322_r2)^ system. In this case, the magnetic interaction is sensitive to the occupancy 
of the Mn d3j,2_r2 orbital. It is also shown that the magnetic interaction between the interfacial 
Mn layer and the bulk region can be changed. Manganite-based heterostructures thus show a rich 
magnetic behavior. We also present how to generalize the argument including t2g orbitals. 
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Transition metal (TM) oxides have been one of the 
main subjects of materials science for decades. Experi- 
mental and theoretical efforts are driven by their rich, 
complex, and potentially useful behaviors originating 
from strong correlations between electrons and/or elec- 
trons and latticesii The recent developments in the crys- 
tal growth techniques, in particular the (laser) molecular- 
beam epitaxy, have made us recognize the opportunity to 
further control their behaviors and to generate phenom- 
ena that are not realized in the bulk systems i-^-^- 

Here, we focus on the magnetic behavior at an inter- 
face between perovskite manganite and other TM ox- 
ides. Perovskite manganites, especially Lai_2;Sr2:Mn03 
(LSMO), are particularly important because of their fer- 
romagnetic (F) metallic behavior with relatively high 
Curie temperature Tc and large polarization. Control- 
ling the magnetic interaction at interfaces involving man- 
ganites would cause a technological breakthrough for 
electronic devices using, for example, a tunneling magne- 
toresistance (TMR) effecli^'^ and an exchange bias (EB) 
effect .-iS.This requires the microscopic information on the 
orbital states, not only on the spin states as demonstrated 
for cuprate/manganite interfaces in Refs.lolandllJI. How- 
ever, it remains controversial whether the magnetic mo- 
ment is induced in the cuprate regiort^iil or the dead 
layers appear in the manganite region^ii^ 

The difficulty dealing with interfaces involving 
strongly-correlated electron systems comes from the 
small volume fraction which makes the experimental 
analysis challenging, and strong-correlation effects which 
hinder some of theoretical treatments. Therefore, if a 
Goodcnough-Kanamori-typei^^ transparent description 
of the interfacial magnetic interaction becomes available, 
both experiment and theory would greatly benefit. 

In this paper, we present a general consideration for the 
magnetic interaction at an interface involving mangan- 
ites. We first focus on the interfacial interaction derived 
by d3z2_r^ orbitals which have the largest hybridization 



along the z layer-stacking direction. We see that the sign 
of the magnetic interaction via the ^322^^2 orbitals is nat- 
urally fixed based on the molecular orbitals formed at the 
interface and the generalized Hund's rule. The argument 
uses localized orbitals, and therefore shows only the qual- 
itative trend. The molecular orbitals effectively lift the 
degeneracy between d^y2_j.2 and dx^~y^ orbitals by the 
order of the hopping intensity. In the second part, we per- 
form the model Hartree-Fock calculation and show that 
the broken degeneracy can lead to the additional change 
in the magnetic interaction between the interfacial Mn 
layer and its neighboring Mn layer. We also discuss how 
to generalize the molecular-orbital based argument for 
more complicated situations including t2g orbitals. 



II. MOLECULAR-ORBITAL PICTURE 

In this section, we consider the magnetic interaction 
between manganite and other TM oxides focusing on the 
molecular orbitals formed by d^z^_r2 orbitals which have 
the largest overlap at the interface. The TM region is 
specified by the number of electrons occupying a d3z^_j.2 
orbital. Here, t2g electrons are assumed to be electroni- 
cally inactive and considered as localized spins when fi- 
nite number of electrons occupy i2g orbitals. Generaliza- 
tion including these electrons will be discussed later. 

(•^Sz^-r^)" system. Let us start from the simplest case, 
an interface between Mn and a {d3z2_r2)'^ system (Fig.[T|). 
In this case, the bonding (B) orbital is occupied by an 
electron whose spin is parallel to the localized i2g spin 
in Mn, while the antibonding (AB) orbital is unoccu- 
pied. When there are other unpaired electrons in the 
(d^z'^ _j.2)'^ system at dx2-y2 and/or i2g orbitals, their 
spins align parallel to that of the electron in the B or- 
bital due to the Hund coupling. Thus, the F coupling is 
generated between Mn and (d^z^-r^)'^ systems. This is 
equivalent to the double-exchange (DE) interaction orig- 
inally proposed by Zener.-'^^ When d^2_y2 is much lower 
in energy than d^z^_r2 and in the interfacial Mn layer 
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FIG. 1: (Color online) Molecular orbital (middle) formed by 
d3^2_^2 orbitals on Mn (right) and TM with the {d^^2_^2) 
configuration (left). 



(termed d^2_y2 order), Mn and (d^z'^-r^)^ systems are 
virtually decoupled. Thus, the magnetic coupling is due 
to the superexchange (SE) interaction between t2g elec- 
trons. This interaction is either F or antiferromagnetic 
(AF) depending on the orbital state and the occupancy 
of the TM t2g level. 

(d^z^-r^)^''^ normal. This simple consideration can be 
easily generalized to {d^z^-r^y and (d^z'^-r^)'^ systems. 
First we consider that the d^z^-r^ and d^^-y^ are nearly 
degenerate in the interfacial Mn layer and the unoccu- 
pied d^2_y2 level in the TM region is much higher than 
the djz^-r^ level (Fig. [51 top figures). We call this config- 
uration "normal" (N) configuration. In the lowest energy 
configuration, B orbitals and the Mn d^2_y2 orbital are 
occupied by electrons. For the {d^z'^-r^Y system, the F 
interaction is favorable as in the (d^z^-r^y^ system. On 
the other hand, for the {d^^z^-r'^Y system, the down elec- 
tron orbital is hybridized with the minority band in the 
Mn region. Thus, the "down" B orbital is higher in en- 
ergy and has larger weight on the TM than the "up" B 
orbital. Because of the Hund coupling with the "down" 
electron in the B orbital, other unpaired electrons, if they 
exist in d,j.2_y2 and/or t2g orbitals, tend to be antiparallel 
to the Mn spin. 

Since dx2_y2 orbitals are predominantly occupied in 
the interfacial Mn layer due to the B/AB splitting of 
(i3z2_r2-based molecular orbitals, further stabilization of 
dr^2_y2 orbitals in the Mn region, i.e., d.j.2_y2 order, would 
not affect the interfacial magnetic coupling discussed 
here. But, this could reverse the magnetic coupling be- 
tween the interfacial Mn layer and the second Mn layer 
as discussed in the next section. 

id^z'^-r'^Y''^ anomalous. When the d.j;2_y2 level in the 
TM region becomes lower than the Mn dj.2_y2 level, the 
charge transfer occurs. We shall call this configuration 
"anomalous" (AN) configuration (Fig. [31 middle figures). 
The electron transferred to the TM d^'^-y'^ orbital has the 
same spin as the higher energy B orbital due to the Hund 
coupling (indicated by arrows). Therefore, the sign of 
the magnetic coupling between the Mn and {d^z^-r'^Y''^ 
systems is unchanged. 

Note that this argument is applicable when the hop- 
ping probability between dx2-y2 orbitals in the Mn and 
the TM regions is negligibly small. The finite hop- 
ping probability would make the charge transfer contin- 
uous. Furthermore, when the hopping probability be- 




FIG. 2: (Color online) Molecular orbitals formed by d3^2_j.2 
orbitals on Mn and the {d^^2_^2y system (left column) and 
the {d^^2_^2)^ system (right column). In the normal (anoma- 
lous) configurations, d3^2_^2 and d^2_y2 orbitals are nearly 
degenerate in the interfacial Mn, and the unoccupied d^2_y2 
orbital in the neighboring TM is higher in energy (lower in 
energy than the occupied Mn d^2_y2). In the JT case, the 
d322„^2 level is significantly lower than the d^2_y2 level. Black 
(light) lines indicate the level of majority (minority) spins. 
The up level and down level are exchange split resulting in 
the level scheme as indicated. The minority levels are ne- 
glected in the upper left two because these are irrelevant. 



comes large, the DE interaction is generated. Although 
the DE interaction through d.j.2_y2 bonds may not be re- 
alistic, it cooperatively stabilizes the F spin alignment 
for the {d^z^-r^Y case, while it competes with the AF 
tendency for the {d^z^_^2)'^ case. 

Magnetic interactions discussed so far are insensitive 
to the electron density in the interfacial Mn because the 
interactions are mainly derived from the virtual electron 
excitation from the occupied d^z^-r^ orbital in the TM 
region to the unoccupied counterpart in the Mn region. 

Next, we consider that the d^z^-r^ level is much lower 
than the dx2-y2 in the interfacial Mn layer due to ei- 
ther the local Jahn- Teller distortion or compressive strain 
originating from the substrate (Fig. [51 lower figures de- 
noted by "JT"). The magnetic coupling in this case is 
sensitive to the electron density of the interfacial Mn. 

(d3^2_^2)^ JT. When the Mn d^z^-r^ density is close 
to 1, the SE interaction between the occupied d^z^-r^ or- 
bitals becomes AF. On the other hand, when the density 
is much less than 1, the F interaction between {d^z'^-r^Y 
configuration on Mn and (d^z^-r^Y becomes dominant. 

{d^z^_^2)'^ JT. When the Mn d^z'^-r^ occupancy is 
close to 1, "up" electrons are localized on each sites be- 
cause both B and AB molecular orbitals are occupied 
while down electrons can be excited or leaked from the 
(d^z'^-r^Y system to the Mn minority level, i.e., down 
electron density is virtually reduced in the TM region. As 
a result, unpaired spins, if they exist in dj.2_y2 and/or t2g 



orbitals, become parallel to the up spin, i.e., F coupling. 
When the Mn d^z2_r2 density becomes much less than 1, 
the up AB orbital becomes less occupied while keeping 
the occupancy of B orbitals relatively unchanged. Even- 
tually, the down density in the TM region becomes larger 
than the up density, and the magnetic coupling between 
the Mn and {d^z2_j.2)'^ regions becomes AF. 



III. MODEL HARTREE-FOCK ANALYSIS 

In the previous section, we discussed the interfacial 
magnetic coupling controlled by the molecular orbitals. 
The separation between B and AB molecular levels can 
become as large as the order of t, the hybridization be- 
tween fi3z2_^2 orbitals along the z direction. Since the 
interfacial Mn d:^z2_r2 band is represented by the B (AB) 
(i322_^2 orbital for a (d322_^2)''(-'^'^'/manganite interface, 
the Cg degeneracy is effectively lifted in the interface 
layer. This degeneracy lifting is expected to affect the 
magnetic interaction in the Mn region. In this section, 
we discuss this effect using the microscopic model calcu- 
lation. 

We consider a two-band DE model given by 

H = Y, A"»a - J2 {^'^U^Aadjb + h.c.) 

i {ij)ab 

+ ^Unianii3 + j'^Sti- Stj. (1) 



Here, an electron annihilation operator at site i and or- 
bital a[= a{d322_r2),/3{dx2-y2)] is given by dia, n,a = 
d\^dia, and the level difference between a and /3 is given 
by A. We consider the large Hund coupling limit, in 
which the spin direction of a conduction electron is al- 
ways parallel to that of a localized t2g spin on the same 
site, and omit the spin index. Instead, the relative orien- 
tation of t2g spins is reflected in the hopping matrix; Uij 
is the unitary transformation representing the rotation 
of the spin direction between sites i and j . For simplic- 
ity, we only consider nearest-neighboring (NN) hoppings 
between Mn eg orbitals via oxygen 2p in the middle. Us- 
ing the Slater-Koster scheme, ^^ the orbital dependence 
of f^ is written as i^f+, = ^;^^^^^ = t, if^f^^^^^ = f , 

^ _ +/3" _ (_\VM and /"^■'5" - fPP _ n 

i,i+x(y) ~ ^i,i+x{y) — \ I i: ^ ^^^" ^i,i+z — ''i.i+z — "' 

The third term represents interorbital Coulomb inter- 
action. Due to the eg symmetry, U is related to the 
intraorbital Coulomb interaction U and the interorbital 
exchange integral Jjf as U — U — 3Jh- The last term 
represents the AF SE interaction between NN t2g spins 

\St\ = l 

From the optical measurements, the on-site inter- 
actions are estimated as C/ ^ 3 eV and Jh ^ 
0.5 eV.— The density-functional theory calculation pro- 
vides t ~ 0.5 eV<^ Using the mean-field analysis for 
the Neel temperature ^ 120K of CaMnOa, one esti- 
mates J '^ 1 meV.— A similar value is obtained from the 
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magnon excitations in the A-AF phases of 50 % doped 
Pri_2:Sr2;Mn03 and Ndi-^jSr^^MnOa supposing that the 
AF interaction is due to the same J^ Thus, in what fol- 
lows, we take U = 3i. Considering some ambiguity, the 
realistic value for JSf/t is expected to be r^ 0.01-0.05. 

We analyze the model Hamiltonian, Eq. ([1]), using the 
Hartree-Fock approximation at T = focusing on the 
doped region (carrier density N far away from 1). In light 
of the experimental reports, we compare the energy of 
the following eight magnetic orderings: F ordering, pla- 
nar AF ordering in which spins align ferromagnetically in 
the xy {xz or yz) plane [A (A')], chain-type AF ordering 
in which spins align ferromagnetically along the z (x or 
y) direction [C (C)], zigzag AF in which spins form fer- 
romagnetic zigzag chains in the xy {xz or yz) plane [CE 
(CE')], and NaCl-type AF (G). At iV -^ 1, in addition 
to the spin symmetry breaking, orbital symmetry can be 
broken due to the SE mechanism in the present model;^ 
Since we are focusing on the metallic regime A'^ < 1, we 
do not consider such a symmetry breaking. 

The numerical results for the bulk phase diagram are 
presented in Figs.|3](a)-[3](c). Here, all phase boundaries 
are of first order, and those at small A^ can be replaced 
by canted AF phases or the phase separation between the 
undoped G-AF phase and doped F or AF phases. The 
overall feature is consistent with the previous theoretical 
reports..2^--22. At A = 0, A-AF and A'-AF (C- and C'-, 
CE- and CE'-) are degenerate but the degeneracy is lifted 
by the finite A. We found the CE phase at C/ = A = at 
JSl ^ 0.112i and N > 0.5 (not shown) as in the previous 
reportsi^"— At {/ = 3i, the CE phase becomes unstable 
against A- and C-AF phases and appears only at the 
positive A with JS^ ^ O.li. JS'^ ~ 0.05i reproduces the 
phase diagram of a high Tc system such as LSMO and 
Pri_3;Sr3;Mn03 (Refs. IH and Isl fairly well. 

The main effect of the level separation A is chang- 
ing the stability of planar-type AF (A or A') with re- 
spect to the chain-type AF (C or C) and F states. In 
particular, the A-AF phase is stabilized by the positive 
A more strongly than the C-AF phase by the negative 
A. This is because the energy gain by the DE mech- 
anism is favorable for the A-AF than the C-AF. The 
result is semiquantitatively consistent with the previ- 
ous report based on the density-functional theory<^ At 
A = t, the boundary between F and A-AF phases is 
moved down to JS"^ ~ 0.02t at 0.3 < A^ < 0.7. This be- 
havior suggests that, when the dy,z2-r^ AB level for the 
((i322_r2)^'^/manganite interfaces is about t higher than 
the Mn d^2_y2 level, the magnetic coupling between the 
interfacial Mn and the second Mn layers is switched to 
AF while retaining the intraplane F coupling. 

We confirmed this behavior by computing the surface 
phase diagram considering F phase and two AF phases: 
Al(2) where the surface layer (and the second layer) is 
antiferromagnetically coupled to its neighbor. We intro- 
duce positive A only on the surface layers in the 20-layer 
slab. As shown in Fig. 131(d), a large part of F phase 
is replaced by Al phase compared with the bulk phase 




FIG. 3: [(a)-(c)] Mean-field phase diagrams of doped mangan- 
ites as a function of electron density A^ and the AF interaction 
J for three choices of the level difference between eg orbitals 
A. At A < (>)0, d3z2_r2 is lower (higher) in energy than 
d^2_y2 . For notations of the magnetic phases, see the main 
text, (d) Phase diagram for the 20-layer slab with A = i in 
the surface layers. A*' in this case corresponds to the mean 
electron density. For Al and A2 phases, schematic spin align- 
ments are also shown. A dashed line is the phase boundary 
between F and A-AF phases in the bulk calculation. 



diagram (b). (Precise phase boundary requires detailed 
information of the surface or interface.) Although the 
parameter regime is small, it is also possible that surface 
three layers are AF coupled while the other couplings re- 
main F, A2 phase, before the whole system enters A-AF 
when J is increased or N is decreased. When the d:}z^_r^ 
orbital is stabilized, the inter-Mn-layer coupling remains 
F but the intraplane F coupling is reduced. Therefore, 
in-plane canted AF structure may result for small N. 



IV. SUMMARY AND DISCUSSION 

Summarizing, we presented a general consideration on 
the magnetic interaction between the doped mangan- 
ite and other transition metal oxides when an interface 
is formed. Using the molecular orbital formed at the 
interface and the generalized Hund's rule, the sign of 
the magnetic interaction is determined (Sec. II). The 
bonding/antibonding splitting of the molecular orbitals 
leads to the degeneracy lifting of Cg orbitals on the in- 
terface Mn layer. Further, the bulk strain lifts the eg 
degeneracy. These effects control the magnetic interac- 
tion in the interfacial Mn plane and between the inter- 
facial Mn plane and its neighbor (Sec. III). Consider- 



TABLE L Magnetic interaction at an interface between Mn 
and TM with the {d^^2_^2)" configuration. The interfacial 
Mn is indicated by Mn(l), and Mn in the second layer by 
Mn(2). Mn(l)-Mn(l) indicates intraplane interaction while 
the others interplane interactions. The TM-Mn(l) interac- 
tion is based on the molecular-orbital picture presented in 
Sec. II while the Mn-Mn interaction is based on the model 
Hartree-Fock study presented in Sec. III. At d^2_y2 order, 
d^2_y2 orbital is stabilized at the (interfacial) Mn layer. At 
F*, Mn-Mn interaction is weak, and the canted AF ordering 
may result. See the stabilization of the C-AF phase by the 
JT-type distortion A < in Fig. [3] (a), the stabilization of 
the C- and A-AF phases by reducing the carrier density A*' 
in Fig. |3] (b), and the stabilization of the Al phase by the 
interfacial d^2_y2 order in Fig. |3](d). 



n Condition TM-Mn(l) Mn(l)-Mn(l) Mn(l)-Mn(2) 
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F* 
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AF 
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AF 
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AF 
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AF 
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F 


F* 
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JT w/ small A^ 


AF 


F* 
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ing these effects, we summarized the magnetic couplings 
in (d322_r2)"/manganite interfaces in TableU) Although 
the present argument is rather qualitative, it is physically 
transparent and can be applied to a variety of systems. 
It is also straightforward to generalize the argument to 
include other orbitals. Therefore, the present argument 
will also help a more quantitative analysis with detailed 
information from either the experiment or the first prin- 
ciple theory. 

It is worth discussing the implication of the present re- 
sults to the real systems. An example of the {d3z2_r2)'^ 
system is high- Tc cuprate. It has been reported that 
the magnetic coupling between YBa2Cu307 (YBCO) and 
Lai_a:Caa;Mn03 is AF)i and d3^2_j.2 and dx2_y2 in the 
interfacial Cu have a similar amount of holes.— This cor- 
responds to the AN situation. F coupling due to the 
DE remains in the Mn region because of the finite band 
width of dj.2_y2. An example of the {d^z^_j.2)^ system 
is BiFeOa (BFO). Recently, the EB effect was reported 
at BFO/LSMO interfaces accompanying the "AF" cou- 
pling between BFO and LSMO.^^ We expect the N sit- 
uation with dj.2_y2 ordering at this interface. Although 
the interfacial coupling is F, the AF coupling between 
the interfacial Mn and the second Mn layers results in 
the AF alignment between BFO and bulk LSMO as ob- 
served experimentally and is responsible for the exchange 
bias effect. 



A question one may ask is what causes the "AN situ- 
ation" in YBCO and the "N situation" in BFO? A qual- 
itative explanation is as follow: in YBCO, the unoccu- 
pied Cu d^2_y2 state is right above the Mott gap and 
its position is nearly identical to the occupied band of 
manganites.— Therefore, the charge transfer from Mn eg 
to Cu dx2-y2 can easily occur. On the other hand, the 
high-spin state is realized in BFO, and the unoccupied 
dx2-y2 state with opposite spin with respect to the ma- 
jority electrons is located far above the gap. In addition, 
the very close chemical potentials (i.e. close d^^2_^2 lev- 
els) of BFO (Ref. Hi) and LSMO (Ref. HI) maximize the 
B and AB splittings. This situation is favorable for the 
d^2_y2 ordering in the interfacial Mn layer and the re- 
sulting AF coupling between the first and the second Mn 
layers [see Fig. |3] (d)] . 

Finally, an example of the (d^z^^^^Y' system may 
be non- magnetic SrTiOa, and the coupling with this 
is expected to affect the magnetic state near the in- 
terfacial Mn. For small doping x of LSMO, the cou- 
pling with SrTiOs (with the smaller lattice constant of 
SrTiOs) increases the d^z'^-r'^ orbital occupancy sup- 
pressing the inplane DE effect. For large doping, SrTiOs 
creates the tensile strain stabilizing dx2_y2, and the out- 
of-plane F coupling is reduced. ^•^ Both are expected to 
cause a more rapid decrease in the ordered moment with 
increasing temperature than in the bulk region^ re- 
sulting in the rapid suppression of the TMR effect in 
LSMO/SrTiOs/LSMO junctions J^^^ The in-plane (out- 
of-plane) spin canting may also be realized in the for- 
mer (latter). ^^ For undoped LaMnOa, the out-of-plane 
ferromagnetic coupling may result because the overlap 
between the occupied d3j,2_^2 in the first Mn layer and 
the unoccupied dx2^z^ or dy2_z^ orbitals in the second 
Mn layer is increased, favorable for the F SE interaction 
between Mn layers. 

However, since t2g orbitals in titanates are located near 
(slightly above) the Fermi level of manganite,— one may 
need to consider t2g orbitals more carefully as discussed 
below. 

Extention to t2g systems. In t2g systems such as ti- 
tanates, vanadates, and cromates, coupling between t2g 
orbitals could become as important as the coupling be- 
tween eg orbitals. Here, we discuss how to generalize 
the molecular-orbital argument presented in Sec. II to 
t2g systems. As an example, we consider an interface be- 
tween titanate with the d*^ configuration and manganites. 
Extending the argument to other systems is straightfor- 
ward. 

Figure S] (a) shows the level diagram of the ti- 
tanate/manganite interface including both eg and t2g 
orbitals. For simplicity, only bonding orbitals are pre- 
sented. Because the d3z^_r^ level in titanate is far above 
the occupied levels in manganite and the unoccupied t2g 
levels in titanate and manganite (down electrons for the 
latter) are close?^ highest occupied molecular orbitals 
could be either B up d^^^^-r^ or d^2_y orbital [Fig. |3](b) 
which is equivalent to Fig. [Tj or B down dxz,yz orbitals 
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FIG. 4: (Color online) Molecular orbitals formed by 3d or- 
bitals on Mn and Ti, originally dP . Here, only bonding or- 
bitals are shown, (a) Full level diagram including both eg 
and t2g- Black (light) lines indicate the level of majority (mi- 
nority) spins. The highest occupied molecular orbitals and 
the magnetic alignment depend sensitively on the detail of 
the interface as shown in (b) and (c). (b) [(c)] The up B or- 
bital of d^^2_^2 orbital is lower (higher) in energy than the 
down B orbitals of dxz,yz- In (b), induced magnetic moment 
in Ti is parallel to Mn, i.e., {d^^2_^2)'^ configuration, while in 
(c), it depends on the relative occupancy of Ti dxz,yz orbitals 
in the up and down B orbitals. When the occupancy of the 
Ti down dxz,yz orbitals is larger than the up d^z.yz orbitals, 
net moment induced in Ti site becomes antiparallel to the Mn 
moment. 



[Fig. 2] (c)]. Note that the interfacial hybridization be- 
tween dxy orbitals on Ti and Mn is much smaller than 
those between d^z and between dyz- 

In the case of Fig. [4] (b), induced moment in the ti- 
tanate region is tiny but parallel to the moment in man- 
ganite region. On the other hand in the case of Fig. 0] 
(c), the induced moment in the titanate region could be 
either parallel or antiparallel to the manganite moment. 
This depends on the relative weight of the down electron 
density in the B dxz,yz orbitals with respect to the up 
electron density in the B dxz,yz orbitals. The situation 
in Fig. |4] (c) with antiparallel spin arrangement between 
titanate and manganite could happen when the original 
electron density in the manganite eg orbital is large and 
d^z^-r^ level both in titanate and manganite is high due, 
for example, to the in-plane tensile strain. But in gen- 
eral, the difference between two configurations. Figs. 0] 
(b) and 2] (c) with either parallel or antiparallel spin 
configurations, would be subtle. Therefore, depending 
on a variety of condition such as the sample preparation, 
any situation could be realized. 

When the number of t2g electrons is increased, such as 
doped titanates, vanadates, and cromates, electrons tend 
to enter the down B dxz,yz orbitals and, then, the down 
B dxy, resulting in the antiparallel spin configuration. 
However, the antiparallel configuration becomes unsta- 
ble against the parallel configuration when the electron 



number in t2g orbitals becomes large and the level sep- 
aration between t2g and Cg orbitals, i.e., \QDq, becomes 
relatively small. 

In this case, because of the strong on-site Coulomb 
interactions, the energy gain by forming B orbitals be- 
comes small for t2g electrons and comparable to having 
electrons in both B and AB orbitals with the parallel spin 
configuration. The parallel configuration further lowers 
the energy by forming d^z^^r^ B orbital and by the Hund 
coupling between the d3z2_j.2 B orbital and t2g electrons 
on the t2g system, i.e., the Zener's double-exchange fer- 
romagnetism discussed in Sec. II. 

When there are more than three d electrons with rel- 
atively large lODq, a low spin state is realized. In this 
case, three electrons enter B orbitals formed with t2g mi- 
nority bands of Mn and remaining electrons enter AB 
orbitals formed with f2g majority band of Mn. Thus, the 
antiparallel configurations persist. Such a situation may 
be realized in, for example, an interface between man- 
ganites and SrRuOs in which Ru'*"'" is in a low spin state 
with t2g- This AF configuration can be turned to the 
F configuration when the double-exchange-type interac- 
tion becomes dominant due to the formation of d3z2_r2 
B orbitalJI 



So far, we have considered the ideal lattice structure in 
which orbitals with different symmetry do not hybridize. 
In reality, bond angle formed by two transition metal ions 
and an oxygen ion in between becomes smaller than 180° 
allowing electrons to hop between orbitals with different 
symmetry. As a result, additional magnetic channels are 
generated. A simple argument presented in this paper 
can be generalized to deal with such a situation. 



Acknowledgments 

The author thanks P. Yu, R. Ramesh, J. Santamaria, 
and C. Panagopoulos for stimulating discussions and 
sharing the experimental data prior to publication, J. 
Kunes for discussion, and the Kavli Institute of Theo- 
retical Physics, University of California Santa Barbara, 
which is supported in part by the National Science Foun- 
dation under Grant No. PHY05-51164, for hospitality. 
This work was supported by the Materials Sciences and 
Engineering Division, Office of Basic Energy Sciences, 
U.S. Department of Energy. 



^ M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 

70, 1039 (1998). 
^ M. Izumi, Y. Ogimoto, Y. Konishi, T. Manako, M. 

Kawasaki, and Y. Tokura, Mater. Sci. Eng., B 84, 53 

(2001). 
^ A. Ohtomo, D. A. Muller, J. L. Grazul, and H. Y. Hwang, 

Nature (London) 419, 378 (2002). 

* S. Okamoto and A. J. Millis, Nature (London) 428, 630 
(2004). 

^ A. Ohtomo and H. Y. Hwang, Nature (London) 427, 423 
(2004); A. Brinkman, M. Huijben, M. van Zaik, J. Huijben, 
U. Zeitler, J. C. Mann, W. G. van der Wiel, G. Rijnders, 
D. H. A. Blank, and H. Hilgenkamp, Nature Mater. 6, 
493 (2007); N. Reyren, S. Thiel, A. D. Caviglia, L. Fit- 
ting Kourkoutis, G. Hammerl, C. Richter, C. W. Schnei- 
der, T. Kopp, A.-S. Ruetschi, D. Jaccard, M. Gabay, D. 
A. Muller, J.-M. Triscone, and J. Mannhart, Science 317, 
1196 (2007). 

^ J. Stahn, J. Chakhalian, C. Niedermayer, J. Hoppler, T. 
Gutberlet, J. Voigt, F. Treubel, H.-U. Habermeier, G. 
Cristiani, B. Keimer, and C. Bernhard, Phys. Rev. B 71, 
140509(R) (2005). 

'^ J. Chakhalian, J. W. Freeland, G. Srajer, J. Strempfer, 
G. Khahullin, J. C. Cezar, T. Charlton, R. Dalgliesh, C. 
Bernhard, G. Cristian, H.-U. Habermeier, and B. Keimer, 
Nat. Phys. 2, 244 (2006). 

* A. Hoffmann, S. G. E. te Velthuis, Z. Sefrioui, J. Santa- 
maria, M. R. Fitzsimmons, S. Park, and M. Varela, Phys. 
Rev. B 72, 140407 (R) (2005). 

® J. Chakhalian, J. W. Freeland, H.-U. Habermeier, G. Cris- 
tiani, G. Khaliullin, M. van Veenendaal, and B. Keimer, 
Science 318, 1114 (2007). 
"^ A. Bhattacharya, S. J. May, S. G. E. te Velthuis, M. Waru- 



sawithana, X. Zhai, B. Jiang, J.-M. Zuo, M. R. Fitzsim- 
mons, S. D. Bader, and J. N. Eckstein, Phys. Rev. Lett. 

100, 257203 (2008). 
" M. Bowen, M. Bibes, A. Barthelemy, J. P. Contour, A. 

Anane, Y. Lemaitrem, and A. Fert, Appl. Phys. Lett. 82, 

233 (2003). 
^^ Y. Ogimoto, M. Izumi, A. Sawa, T. Manako, H. Sato, H. 

Akoh, M. Kawasaki, and Y. Tokura, Jpn. J. Appl. Phys., 

Part 2 42, L369 (2003). 
" P. Yu, J.-S. Lee, S. Okamoto, M. D. Rossell, M. Huijben, 

C.-H. Yang, Q. He, J.-X. Zhang, S. Y. Yang, M. J. Lee, 

Q. M. Ramasse, R. Erni, Y.-H. Chu, D. A. Arena, C.-C. 

Kao, L. W. Martin, and R. Ramesh, Phys. Rev. Lett. 105, 

027201 (2010). 
^* M. van Veenendaal, Phys. Rev. B 78, 165415 (2008). 
^^ W. Luo, S. J. Pennycook, and S. T. Pantelides, Phys. Rev. 

Lett. 101, 247204 (2008). 
^® J. B. Goodenough, Magnetism and the Chemical Bond (In- 

terscience, New York, 1963). 
" J. Kanamori, J. Phys. Chem. Solids 10, 87 (1959). 
^® C. Zener, Phys. Rev. 82, 403 (1951). 

^^ J. C. Slater and G. F. Koster, Phys. Rev. 94, 1498 (1954). 
^° N. N. Kovaleva, A. V. Boris, C. Bernhard, A. Kulakov, A. 

Pimenov, A. M. Balbashov, G. Khaliullin, and B. Keimer, 

Phys. Rev. Lett. 93, 147204 (2004). 
^^ C. Ederer, C. Lin, and A. J. Millis, Phys. Rev. B 76, 

155105 (2007). 
^^ E. O. WoUan and W. C. Koehler, Phys. Rev. 100, 545 

(1955). 
^^ H. Kawano-Furukawa, R. Kajimoto, H. Yoshizawa, Y. 

Tomioka, H. Kuwahara, and Y. Tokura, Phys. Rev. B 67, 

174422 (2003). 
2'' K. I. Kugel' and D. I. Khomskii, Sov. Phys. Usp. 25, 231 



(1982). 
^ R. Maezono, S. Ishihara, and N. Nagaosa, Phys. Rev. B 

58, 11 583 (1998). 
^^ J. van den Brink and D. Khomskii, Phys. Rev. Lett. 82, 

1016 (1999). 
^"^ Z. Fang, I. V. Solovyev, and K. Terakura, Phys. Rev. Lett. 

84, 3169 (2000). 
^* L V. Solovyev and K. Terakura, Phys. Rev. Lett. 83, 2825 

(1999). 
^^ J. van den Brink, G. Khahulhn, and D. Khomskii, Phys. 

Rev. Lett. 83, 5118 (1999). 
^° L. Brey, Phys. Rev. B, 71, 174426 (2005). 
^^ Z. Jirak, J. Hejtmanek, E. Pollert, C. Martin, A. Maignan, 

B. Raveau, and M. M. Savosta, J. Appl. Phys. 89, 7404 

(2001). 



^^ O. Chmaissem, B. Dabrowski, S. Kolesnik, J. Mais, J. D. 
Jorgensen, and S. Short, Phys. Rev. B 67, 094431 (2003). 

^^ S. Yunoki, A. Moreo, E. Dagotto, S. Okamoto, S. S. Kan- 
charla, and A. Fujimori, Phys. Rev. B 76, 064532 (2007). 

^'' P.-F. Paradis, T. Ishikawa, and S. Yoda, Meas. Sci. Tech- 
nol. 16, 452 (2005). 

^5 M. P. de Jong, V. A. Dediu, C. TaUani, and W. R. Sala- 
neck, J. Appl. Phys. 94, 7292 (2003). 

^^ S. Okamoto, J. Phys.: Condens. Matter 21, 355601 (2009). 

^■^ J. W. Seo, W. Prellier, P. Padhan, P. BouUay, J.-Y. Kim, 
H. G. Lee, C. D. Batista, I. Martin, E. E. M. Chia, T. Wu, 
B.-G. Cho, and C. Panagopoulos, arXiv: 1006.3603 (unpub- 
lished). 



